Interlayer superfluidity in bilayer systems of fermionic polar molecules 
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We consider fermionic polar molecules in a bilayer geometry where they are oriented perpendicu- 
larly to the layers, which permits both low inelastic losses and superfluid pairing. The dipole-dipole 
interaction between molecules of different layers leads to the emergence of interlayer superfluids. 
The superfluid regimes range from BCS-like fermionic superfluidity with a high T c to Bose-Einstein 
(quasi-)condensation of interlayer dimers, thus exhibiting a peculiar BCS-BEC crossover. We show 
that one can cover the entire crossover regime under current experimental conditions. 
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Ultracold gases of dipolar particles attract great in- 
terest because the dipole-dipole interaction drastically 
changes the nature of quantum degenerate regimes com- 
pared to ordinary short-range interacting gases [l|, Q. 
This has been demonstrated in experiments with Bosc- 
condensed chromium atoms which have a magnetic mo- 
ment of 6/ib equivalent to an electric dipole moment of 
0.05 D [H|. The recent experiments on creatingpolar 
molecules in the ground ro-vibrational state @, @] and 
cooling them towards quantum degeneracy @ have made 
a breakthrough in the field. For such molecules polarized 
by an electric field the dipole-dipole interaction is sev- 
eral orders of magnitude larger than for atomic magnetic 
dipolcs. This opens fascinatingprospects for the obser- 
vation of new quantum phases |ul^.l8l-fl2j. The main ob- 
stacle is the decay of the system due to ultracold chemical 
reactions, such as KRb+KRb=>K2+Rb2 found in JILA 
experiments fl3j | . These reactions are expected to be sup- 
pressed by the intcrmolecular repulsion in 2D geometries 
where the molecules are oriented perpendicularly to the 
plane of their translational motion [14[ . 

In this Letter we consider fermionic polar molecules in 
a bilayer geometry where the dipoles are oriented per- 
pendicularly to the layers (Fig. [J), which leads to low 
inelastic losses and allows for the possibility of superfluid 
pairing. The interaction between dipoles of different lay- 
ers may lead to the emergence of an interlayer superfluid, 
that is a superfluid 2D gas where Cooper pairs are formed 
by fermionic molecules of different layers. We show that 
the interlayer dipole-dipole interaction provides a higher 
superfluid transition temperature than that for 2D spin- 
1/2 fermions with attractive short-range interaction. 

Interestingly, an increase in the interlayer dipole-dipole 
coupling leads to a novel BCS-BEC crossover resembling 
that studied for atomic fermions near a Feshbach reso- 
nance [ill EH ■ The reason is that two dipoles belonging 
to different layers can always form a bound state ljj . As 
long as the binding energy e& is much smaller than the 
Fermi energy Ef, or equivalently the size of the interlayer 



FIG. 1: Bilayer dipolar system under consideration. 

two-body bound state greatly exceeds the intermolecular 
spacing in the {x, y} plane, the ground state of the system 
is the BCS-paired interlayer superfluid. Once a reduction 
of the interlayer spacing A or an increase of the molecu- 
lar dipole moment d by an electric field make tb >> Ep, 
dipolar fermions of different layers form true bound states 
in real space and the ground state is a Bose-condensed 
system of these composite bosons. We describe this pecu- 
liar BCS-BEC crossover and show that interlayer super- 
fluids may be observed for typical parameters of ongoing 
experiments. Strictly speaking, at a finite temperature 
T in the thermodynamic limit this is a crossover from a 
BCS-paired algebraic superfluid to an algebraic bosonic 
superfluid (quasi-BEC) of dimers. However, we keep the 
term BCS-BEC crossover for brevity. 

We consider the bilayer system of Fig. [T] assuming no 
interlayer hopping. The interaction potential between 
two dipoles belonging to different layers has the form: 



V(r) = £ 
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where r is the in-plane separation between these dipolcs. 
The potential V(r) is attractive for r < a/2A, and repul- 
sive at larger distances r. It satisfies the relation 



V(r)d 2 



0. 



(2) 



which precludes the ordinary method of finding a bound 
state in 2D potentials finite at the origin [3]. However, it 



2 



has been proven that V(r) always has a bound state 17| . 
at any dimensionless strength j3 = r*/A, with r* = 
md 2 /h 2 being the dipole-dipole length. For j3 <C 1 the 
binding energy is exponentially small [l9j : 

e b ~ E exp [-8(1 - /3)//3 2 - (5 + 2 7 - 2 In 2)] , (3) 

where i?o = fr 2 /m\ 2 , and 7 = 0.5772 is the Euler con- 
stant. One finds numerically that Eq. ([3]) is valid up to 
P ~ 1. Note that the unusual dependence on the inter- 
action, eft ~ exp(— 8//3 2 ), is a consequence of Eq. ([2]). 

Dipoles of different layers undergo the 2D s-wave scat- 
tering from each other in the interlayer potential V(r). 
We define the off-shell scattering amplitude as 

/(k,k') = (m/h 2 ) [ exp(-tk'r)V(r)Mr)d 2 r, (4) 



where ipk(r) is the true wavefunction of the relative mo- 
tion with momentum k. The potential V(r) shows a slow 
power law decay ~ 1/r 3 at large distances r. Therefore, 
at low relative momenta k <C r^T 1 , A" 1 and k' <C r" 1 , A" 1 
one has two contributions to the scattering amplitude: 
the contribution from short distances and the so-called 
anomalous contribution from distances r ~ 1/k [H| ob- 
tained using a perturbative approach in V(r). The lead- 
ing short-range and anomalous contributions yield the 
following s-wave part of /(k, k'): 



f(k,k') 
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omitting higher order terms. The short-range (logarith- 
mic) contribution is obtained by putting kl = and pro- 
ceeding along the lines of the 2D scattering theory [18| . 
The fc-dcpcndencc of the s-wave part of ipk at distances 
in the interval r*,A <C r <C k" 1 is given by a factor 
[ln(«;/fc) + i7r/2] _1 , where k depends on the behavior 
of V(r) at small r and in the presence of the weakly 
bound state we have k = y/meb/h [Hj]. The anomalous 
term comes from distances where the motion is almost 
free. We then have F^x) = {x 2 /2)F{\/2, 1/2, 2, x 2 ) + 
F(l/2, —1/2, 1, x 2 ), where F is the hypegeometrical func- 
tion, so i*i(l) = 4/7T. This is valid for k' < k, for k' > k 
one should interchange k' and k. A detailed derivation 
of f(k,k'), including fc 2 -terms, will be given elsewhere. 

The anomalous term in Eq. ([5]) corresponds to attrac- 
tion and so does the logarithmic term if k <C k, i.e. if the 
collision energy is much larger than ej,. Thus, both the 
short-range and anomalous contribution may lead to su- 
perfiuid interlayer pairing. Note that for short-range po- 
tentials, like all interatomic potentials decaying as 1/r 6 , 
only the logarithmic term is present in Eq. (|5|). We will 
show that the anomalous scattering drastically influences 
the supcrfluid pairing. The inlayer dipole-dipole interac- 
tion is repulsive and it simply renormalizes the chemical 
potential. This is valid as long as the inlayer repulsion is 
sufficiently weak to exclude crystallization 1(J U | ■ 



Treating molecules of the first and second layers as 
spin-up and spin-down fcrmions our problem is mapped 
onto spin- 1/2 fermions with a peculiar interaction poten- 
tial. For a weak interlayer attractive interaction we use 
the BCS approach and obtain the standard gap equation 
for the momentum-space order parameter: 



A(k) 



dk' 2 V(k-k')A(k') 
(2^r)2 2^ 



tanh (H), (6) 



where e k = \J{E^ — ^) + |A(fc)| 2 is the gapped dis- 
persion relation, [i is the chemical potential, and Ek = 
h 2 k 2 /2m. We rewrite Eq. © expressing the Fourier com- 
ponent of the interaction potential. VYk — k'), through 
the off-shell scattering amplitude [20(. Assuming that 
the s-wave interaction is the leading channel of supcr- 
fluid pairing the renormalized gap equation reads: 
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f tanh(e fc >/2r) 
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As long as the interaction is really weak and fi ~ 
Ep [ll|, Eq. (0 may be employed for calculating A(fc) 
and the supcrfluid transition temperature. As known, in 
2D the transition from the normal to supcrfluid state is 
of the Kosterlitz-Thouless type. However, in the BCS 
limit the Kosterlitz-Thouless transition temperature is 
very close to the critical temperature T c given by the 
BCS gap equation (22[. Using Eq. (J7J we obtain the re- 
lation between T c and the order parameter on the Fermi 
surface at T = 0, A (fcp), which is the same as in 3D [23j : 



T c = ( e T/7r)A (fcF). 



(8) 



When the short-range logarithmic contribution to 
Eq. ([5]) dominates, the anomalous term can be omitted. 
This is in particular the case for (3 approaching unity and 
sufficiently small values of fci?A. Then, using Eq. we 
recover the well-known results [22 , 24 1 : 



Ao(**0 = V2E F e b ; T c = (e^ /n)y/2E F e b . (9) 

Note that we do not include here the second order 
Gor'kov-Mclik-Barkhudarov corrections. They decrease 
both A (k F ) and T c by a factor of e (25|, but Eq. © 
remains valid. A detailed analysis of the BCS limit up 
to the second order will be given elsewhere. 

For /3 < 1 the anomalous scattering dominates, at 
least for not very low kp. Then the logarithmic term 
in Eq. (|5|) reduces to — 7r/3 2 /2, and it is necessary to 
include quadratic terms in A:. In this case the scatter- 
ing amplitude can be calculated using the second order 
Born approximation. The expression for the off-shell 
amplitude is cumbersome. The on-shell amplitude is 
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f(k, fc)=/'/(l+i/'/4), where f'(k) is real and given by 



f=-8kn 



47r(fcrJ 2 tt/3 2 „ ,,, 
■— ^ tL + 3 ^(fcr,) 2 ln 
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2 
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where the first term is dominant and it follows from the 
second term in Eq. ([5]) at k' = k. 

We now use Eq. to calculate T c . For T -> T c , we 
set A(fc') = in the dispersion relation which becomes 
efe/ = l-Ef — -EfI- The main contribution to the integral 
over dk' comes from the region near the Fermi surface, 
where we put k! = kp in the arguments of A and /, and 
taking k — kp use f(kp, kp) from Eq. (fTO)) . For the rest 
of the integration it is sufficient to use f(kp,k') given 
by the second term of Eq. §5§ and employ the relation 
A(fc) ~ A(k F )f(k, k F )/f(k F , k F ) following from Eq. Q. 
After a straightforward algebra we then find: 



T c = O.lEp — — 



£■0 
E F 



0.46 

exp 



4kpr* 



G{k F \,P)\, (11) 



where G(x,y) = (1 — irx/2 + ny /\&x)~ l . The validity 
of Eq. (fTTj) requires 1 kpX ^> 7T/3/16. One easily 
checks that Eq. (fTTj) gives a significantly higher T c than 
that given by Eq. ([9]) . The numerical solution of Eq. 
also confirms this conclusion for kpX approaching unity. 
Fig. [2] shows that T c strongly departs from Eq. © for 
small (3, an anomalous behavior stemming from the long- 
range character of the interlayer interaction. 

For sufficiently strong interactions [i deviates from Ep, 
and Eq. |7| should be complemented by the number 
equation [26|, [2?} for the 2D density n in one layer: 



d 2 k 

(27T) 2 



Ek - V 



tanh ( ——] 
\2TJ 



(12) 



We found numerically [i and A(fc) from the self-consistent 
solution of Eqs. ([7|) and (fT2|) . Alternatively, we used 
Eq. ((6]) together with (fT2|) . which is adequate since the 
potential V(r) is finite and strongly bounded at small r. 

This approach provides a qualitative description of the 
strongly interacting regime 0, [23] . An increase of e& by 
e.g. increasing leads to bound interlayer dimers when ej, 
becomes much larger than Ep and the chemical potential 
for the fermionic molecules is fi ~ — £b/2, in contrast to 
/i ~ £^i? in the BCS regime. These composite bosons 
condense and we thus have a BCS-BEC crossover. An 
approximate crossover line is marked by the condition 
|i = (Fig. [3]). At sufficiently low T, well above 
this line a dimcr (quasi)BEC occurs whereas well below 
the line the system is a Fermi gas which is superfluid or 
normal, depending on T and density (Fig. [3]). 

For strong interactions, T c calculated from Eqs. ((TJ) 
and (|12p cannot be interpreted as the critical tempera- 
ture for the onset of superfluidity. Instead, it corresponds 
to the temperature of pair dissociation [28j . The temper- 
ature of the Kosterlitz-Thoulcss transition, Tkt, below 
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FIG. 2: Critical temperature T c as a function of the dipole- 
dipole strength j3 for fc^A = 0.5. The numerical solution 
(solid) is higher by at least an order of magnitude than the 
result of Eq. (© (dashed). 



2.4 


Dimer / 




2 - 


BEC/ 










1.6 






1.2 


Superfluid Fermi Gas 










0.8 - 






0.4 - 


Normal Fermi Gas 





0.4 0.6 0.8 



1.2 1.4 



FIG. 3: Phase diagram for T — 0.05Ef, obtained from 
Eqs. and (|12[l . The curves indicate the Kosterlitz- 
Thouless (KT) transition and the n — line. 

which the system is superfluid satisfies the equation [29| : 

T KT = Trh 2 p s (T K T)/2M 2 7 (13) 

where M — 2m is the dimer (Cooper-pair) mass, and 
p s is the superfluid mass density just below Tkt, which 
may be determined from our mean-field equations using 
the known expression for the normal density (23j . In 
Fig. ||wc depict Tkt and T c versus f3 for kpX = 0.5. 
In the BCS regime we retrieve Tkt ~ T c and see that 
the ratio T c /Ef can reach 0.04. For strong interac- 
tions where A(fc^) is a sizable fraction of Ep we ob- 
tain Tkt = 0.125£ F (cf. [30j]). In the intermediate 
regime Tkt interpolates smoothly between these lim- 
its (see Fig. |4]) . On the BEC side of the crossover we take 
into account a noticeable normal fraction, which makes 
T KT lower. We obtain that T KT - 0.1E F for f3 ~ 2.2 and 
very slowly decreases with increasing (3 (decreasing the 
interaction) , in agreement with the result for bosons [3l[ . 
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FIG. 4: Kosterlitz-Thouless transition temperature Tkt and 
the critical BCS temperature T c versus the dipole-dipole 
strength /3 for fcpA = 0.5. 

In conclusion, we have shown that bilayer systems of 
fermionic polar molecules which are expected to have low 
inelastic losses, at the same time may allow the obser- 
vation of interesting regimes of intcrlayer superfluidity. 
These regimes range from fermionic BCS-like superfluid- 
ity with a relatively high T c and Cooper pairs formed by 
molecules of different layers, to quasiBEC of interlayer 
dimers, thus exhibiting a peculiar BCS-BEC crossover. 
For example, by making the interlayer spacing A ~ 250 
nm one achieves fc^A ~ 2 for KRb and LiCs molecules 
at densities n ~ 5 10 8 cm -2 corresponding to Ep ~ 110 
nK. Then, varying the LiCs dipole moment d from 0.35 
to 1.3 D by increasing the electric field to about 1 kV/cm, 
one obtains /3 ranging from 1 to 14 and covers the entire 
crossover regime, with Tjct °f a f° w nanokclvin. For KRb 
molecules the strongly interacting regime can be reached 
for the presently achieved d ~ 0.2 D [l4[ by putting a 
shallow in-plane optical lattice and getting /? > 1 due to 
an increase in the effective mass of molecules. 

Our results open exciting perspectives for future stud- 
ies. Imbalanced Fermi mixtures may be studied by 
preparing layers with different chemical potentials (ef- 
fective magnetic field) or with different densities. An 
increase in the dipole-dipole interaction may lead to 
in-plane Wigncr-like dimer crystallization, and perhaps 
opens routes towards a supersolid dimer gas. 
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Note added: After the completion of this work we 
learned of a recent related work of Potter et al. [32| . 
where interlayer dimerization and superfluidity in mul- 
tistacks of polar Fermi molecules have been considered. 
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